Electronic excitations from a linear-response range-separated hybrid scheme 
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We study linear-response time-dependent density-functional theory (DFT) based on the single- 
determinant range-separated hybrid (RSH) scheme, i.e. combining a long-range Hartree-Fock ex- 
change kernel with a short-range DFT exchange-correlation kernel, for calculating electronic excita- 
tion energies of molecular systems. It is an alternative to the more common long-range correction 
(LC) scheme which combines a long-range Hartree-Fock exchange kernel with a short-range DFT 
exchange kernel and a standard full-range DFT correlation kernel. We discuss the local-density ap- 
proximation (LDA) to the short-range exchange and correlation kernels, and assess the performance 
of the linear-response RSH scheme for singlet — > singlet and singlet — s* triplet valence and Rydberg 
excitations in the N2, CO, H2CO, C2H4, and CeH6 molecules, and for the first charge-transfer ex- 
citation in the C2H4-C2F4 dimer. For these systems, the presence of long-range LDA correlation in 
the ground-state calculation and in the linear-response kernel has only a small impact on the excita- 
tion energies and oscillator strengths, so that the RSH method gives results very similar to the ones 
given by the LC scheme. Like in the LC scheme, the introduction of long-range HF exchange in the 
present method corrects the underestimation of charge-transfer and high-lying Rydberg excitation 
energies obtained with standard (semi)local density- functional approximations, but also leads to 
underestimated excitation energies to low-lying spin-triplet valence states. This latter problem is 
largely cured by the Tamm-Dancoff approximation which leads to a relatively uniform accuracy for 
all excitation energies. This work thus suggests that the present linear-response RSH scheme is a 
reasonable starting approximation for describing electronic excitation energies, even before adding 
an explicit treatment of long-range correlation. 



I. INTRODUCTION 



Range-separated density- functional theory (see, e.g., 
Ref. [0 and references therein) constitutes an alterna- 
tive to standard Kohn-Sham (KS) density-functional the- 
ory (DFT) [2( for ground-state electronic-structure cal- 
culations. It consists in combining wave-function-type 
approximations for long-range electron-electron interac- 
tions with density-functional approximations for short- 
range electron-electron interactions, using a controllable 
range-separation parameter. For example, in the single- 
determinant range-separated hybrid (RSH) scheme Q, 
the long-range Hartree-Fock (HF) exchange energy 
is combined with a short-range exchange-correlation 
density-functional approximation. The long-range cor- 
relation energy is missing in this scheme, but it can be 
added in a second step by many-body perturbation the- 
ory for describing van der Waals dispersion interactions 
for instance [J-Q. A simpler approach is the long-range 
correction (LC) scheme 0, also called RSHX [1(|, which 
consists in applying range separation on exchange only, 
i.e. combining the long-range HF exchange energy with 
a short-range exchange density-functional approximation 
and using a standard full-range correlation density func- 
tional. More complicated decompositions of the exchange 
energy have also been proposed, such as in the CAM- 
B3LYP approximation fill ] - 
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Range separation is also applied in linear-response 
time-dependent density-functional theory (TDDFT) [12( 
for calculating excitation energies and other response 
properties. The first and probably most widely used 
range-separated TDDFT approach is based on the LC 
scheme [13J], and involves a long-range HF exchange ker- 
nel combined with a short-range DFT exchange kernel 
and a standard full-range DFT correlation kernel. It 
has also been proposed to use in this scheme an em- 
pirically modified correlation density functional depend- 
ing on the range-separation parameter [l4[. The CAM- 
B3LYP scheme and other similar schemes have also been 
applied in linear- resp onse theory for calculating excita- 
tion energies [IJ, [L5H26J . In all these schemes, the pres- 
ence of long-range HF exchange greatly improves Ryd- 
berg and charge-transfer excitation energies, in compari- 
son to time-dependent Kohn-Sham (TDKS) calculations 
using standard local or semilocal density-functional ap- 
proximations in which they are strongly underestimated 
(see, e.g., Ref.[27t). 

In this paper, we study a range-separated linear- 
response TDDFT method based on the RSH scheme, i.e. 
combining a long-range HF exchange kernel with a short- 
range DFT exchange-correlation kernel with no long- 
range correlation kernel. The motivation for this range- 
separated TDDFT approach is that, as for exchange, 
the long-range part of standard correlation density- 
functional approximations such as the local-density ap- 
proximation (LDA) is usually inaccurate [l|, [28|, [29| , so 
one may as well remove it. This can be viewed as a first- 
level approximation before adding a more accurate treat- 
ment of long-range correlation, e.g., by linear-response 
density-matrix functional theory (DMFT) [3fJ] or linear- 



response multiconfiguration self-consistent field (MC- 
SCF) theory [3l|. These last approaches are capable of 
describing excited states of double excitation character, 
which are out of reach within a single-determinant linear- 
response scheme using adiabatic exchange-correlation 
kernels (except in a spin- flip formulation [32l. l33l|). 

The main goal of this paper is to test whether the 
range-separated TDDFT method based on the RSH 
scheme is a reasonable starting approximation for cal- 
culating excitation energies of molecular systems, even 
before adding explicit long-range correlations. For this 
purpose, we apply the method to singlet — > singlet and 
singlet — > triplet valence and Rydberg excitations in the 
N 2 , CO, H 2 CO, C2H4, and CeH 6 molecules, and to the 
first charge-transfer (CT) excitation in the C2H4-C2F4 
dimer, and compare with the LC scheme, as well as non- 
range-scparatcd methods. In particular, we study the 
effect of dropping long-range LDA correlation in com- 
parison to the LC scheme. 

The paper is organized as follows. The working equa- 
tions of the linear-response RSH scheme are laid down in 
Section |nl and the short-range DFT exchange and cor- 
relation kernels are discussed in Section HTT1 After giving 
computational details in Section HV[ we report and dis- 
cuss our results in Section [VJ Section Wl\ summarizes our 
conclusions. Technical details are given in Appendices. 
Hartrec atomic units are assumed throughout unless oth- 
erwise indicated. 



II. LINEAR-RESPONSE RANGE-SEPARATED 
HYBRID SCHEME 

A. Ground-state range-separated scheme 

In the RSH scheme [3J, the ground-state energy is 
approximated as the following minimum over single- 
determinant wave functions $, 



E RSH =min{($|T + V ext \<P) + E H [n*] 

(1) 

+ <HF^+CN>™4]}, 



where T is the kinetic energy operator, V ex t is the ex- 
ternal potential operator, -EhM is the Hartree energy 
density functional, 

BhH = 9 / ™( r iM r 2)w ee (|ri -r 2 \)dr 1 dr 2 , (2) 

with the Coulombic electron-electron interaction 
w rk (\ti - r 2 |) = l/|ri - r 2 |, E% HF [$] is the long-range 



HF exchange energy 



^HFW=-^/l( ( I ) K(xi,x 2 )|$)| 2 ^ e (|r 1 -r 2 |)dx 1 dx 2 , 

(3) 
with the one-particle density-matrix operator n 1 (x 1 ,x 2 ) 
and a long-range electron-electron interaction w^Qri — 



r 2 |), and E^ r c [n,m] is the short-range exchange- 
correlation energy functional depending on the total den- 
sity n(r) = nf (r) + nj.(r) and the (collinear) spin mag- 
netization density m(r) = n^(r) — ^4.(1"), written with 
the spin densities n a (r) = n(x) for the space-spin coor- 
dinate x = (r, a) . In this work, the long-range interac- 
tion will be taken as w^(r) = erf (/ir)/r, where the pa- 
rameter [i can be interpreted as the inverse of a smooth 
"cut-off" radius, but other interactions have also been 
considered J34l - l36j . What is neglected in Eq. ((T|) is the 
long-range correlation energy E l J~, but it can be added a 
posteriori by perturbative methods |3j-|8|, [37], |38[ . 

In the LC scheme 0, range separation is applied to 
the exchange energy only and the ground-state energy is 
expressed as 



E LC =min{($|T + V ext \$) + M n t>] 

+ £^ r HF [$] +^ r [n$,m$] +£ ! c[n$,m$]}, 



(4) 



where E* r [n, m] is the short-range exchange energy func- 
tional, and E c [n,m] is the full-range correlation energy 
functional. 



B. Linear-response theory 

Just like in standard TDDFT [12(, time-dependent 
linear-response theory applied to the RSH scheme 
leads to a familiar Dyson-like equation for the 
frequency-dependent 4-point linear response function 
x(xi,x 2 ;x' 1 ,x 2 ;a;) to a time-dependent perturbation 
(dropping the space-spin coordinates for simplicity) 



x- 1 M = xo 1 M-/ff-/^HF-/*" 



(5) 



where Xo(w) is the non- interacting RSH response func- 
tion, fu is the Hartree kernel, 

f H (x 1 ,x 2 ;x' 1 ,x' 2 ) =w ee (|ri -r 2 |)<5(xi - x[)8(x 2 - x 2 ), 

(6) 
fxBF i s the long-range HF exchange kernel, 

/^hf( x i, x 2;xi,x 2 ) = ~w l e r e (\r 1 -r 2 \)S(x 1 ~x' 2 )5(x' 1 -x 2 ), 

(7) 
and /££ is the short-range exchange-correlation kernel 
which is frequency independent in the adiabatic approx- 
imation, 



/^(x 1 ,x 2 ;xi,x 2 ) = /^(xi,x 2 ) ( 5(x 1 -x' 1 )5(x 2 -x 2 ), 



with the 2-point kernel 



/*c(xi,X 2 ) 



S 2 E s x r c [n,m] 
5n(xi)5n(x 2 ) 



(8) 



(9) 



Note that a 4-point formalism is required here because 
of the HF exchange kernel. The excitation energies are 
given by the poles of x( w ) i n u - Working in the basis of 



the RSH spin orbitals {</>fc(x)}, the poles can be found 
by the pseudo-Hermitian eigenvalue problem [39( 



A B 

B* A* 



1 
-1 



(10) 



whose solutions come in pairs: the excitation energy u n 
associated with the eigenvector (X n ,Y n ), and the de- 
excitation energy — u n associated with (Y*,X*). The 
matrix elements of A and B are 



^-iajb — \^a &i)0ij0ab ~r *-*-ia.jbi 



B 



ia,jb 



K 



ia,bj i 



(ii) 



where i, j and a, b refer to occupied and virtual spin or- 
bitals, respectively, Ek is the energy of the spin orbital k, 
and K is the coupling matrix accounting for the contri- 
butions of the different kernels, 



K Mj6 = (aj\f H \ib) + (aj\f* nF \ib) + (aj\f£\ib) 

= (aj\w ee \ib) - (aj\w l e r e \bi) + (aj|/£|»6>, 



(12) 



where (aj\w ee \ib) and (aj \w l e e \bi) are the two-electron in- 
tegrals for the Coulombic and long-range interactions, 
respectively, and {aj\f x r c \ib) are the matrix elements of 
the short-range exchange-correlation kernel, 



(13) 



(ai\fZ\ib) = J C(xi)0*(x 2 ).^(x 1 ,x 2 ) 
x (/>i(xi)0f,(x 2 )dxidx2. 



For real- valued orbitals, and if A + B and A — B are pos- 
itive definite, Eq. (|10[l is conveniently transformed into a 
half-size symmetric eigenvalue equation [39j 



MZ, 



a/2 



B) 



1/2 



(14) 



and the nor- 



-V2 



with M = (A-B) ' (A + B)(A 
malized eigenvectors Z„ = y/uJ^{K — B)~ ±/z (X. + Y n 
The Tamm-Dancoff approximation (TDA) J4Q| consists 
in neglecting the coupling between the excitations and 
the de-excitations, i.e. setting B = 0. We note, in 
passing, that the TDA can also be viewed as a non-self- 
consistent approximation to the static (multiplet-sum) 
ASCF method, which identifies the excited states with 
stationary points on the ground-state e nerg y surface as 
a function of the orbital parameters [3j, EH . 

The same equations apply identically to the LC 
scheme except that the short-range correlation kernel f£ r 
has to be replaced by the full-range one f c \J3^ . 



For simplicity, they will be referred to as "singlet exci- 
tations" and "triplet excitations" . The modifications for 
spin adaptation are located in the expression of the cou- 
pling matrix K, which becomes, for singlet excitations, 



l K ia , jb = 2(aj\w ee \ib) - (aj\w l e r e \bi) + 2(aj\ l f° r c \ib) 



(15) 



and, for triplet excitations, 



3 K iaj5 = -(aj\w l e r e \bi) + 2(aj\ A fZ\ib), (16) 

where the indices i, j, a, b refer now to spatial orbitals and 
the singlet and triplet short-range exchange-correlation 
kernels arc 



! S 2 E s x r c [n,m] 

Jxc\ r li r V - 



and 



7^(ri,r 2 ) 



<5n(ri)<5n(r 2 )' 



(5m(ri)(Sm(r 2 ) ' 



(17) 



(18) 



where the derivatives are taken 
netization density, m(r) = 0. 
dependent exchange functional E* r 
from the spin-independent one E x 
0] via the spin-scaling relation 
(E s J[2n^\ + £'| r [2nJ)/2, one can show that the singlet 
and triplet exchange kernels are identical, and, for closed- 
shell systems, can be written with the spin-independent 
functional, 



at zero spin mag- 
Because the spin- 
n, m] is constructed 
r \n] = E s x r [n,m = 
El, E x r [n,m] = 



/r(ri,r 2 ) = 7r(ri,r 2 ) = 3 f^(r U T a ) = 



FEl 



5n(ri)i5n(r 2 ) 
(19) 
Therefore, contrary to the case of the correlation func- 
tional, the dependence on the spin magnetization density 
does not need to be considered in practice in the exchange 
functional for closed-shell systems. 

The oscillator strength /„ for state n is zero for a triplet 
excitation, and it is calculated with the following formula 
for a singlet excitation (in the dipole length form) [39| 



J2 E^-fC^-'B) 172 ^ 



(20) 
where d a ^ a = J 4'i{^-) r a'Pa(x)dr is the a-component of the 
transition dipole moment between the spatial occupied 
and virtual orbitals </>i(r) and </> a (r). 



C. Spin adaptation for closed-shell systems 

For spin-restricted closed-shell calculations, Eq. (fl4] 
can be decoupled into two independent eigenvalue equa 
tions for singlet —> singlet excitations and for singlet —, 
triplet excitations, respectively |42h44| (see Appcndix [A"| 



III. SHORT-RANGE ADIABATIC 
EXCHANGE-CORRELATION KERNELS 

We will consider here the short-range adiabatic ex- 
change and correlation kernels in the local-density ap- 
proximation (LDA). 



A. Exchange kernel 

The short-range spin-independent LDA exchange en- 
ergy functional is written as 



E 



tc.LDAl 



= er(n(r))dr, 



(21) 



where ej r (n) = ne^, r (n) is the short-range energy den- 
sity defined with the exchange energy per particle e* r (n) 
of the homogeneous electron gas (HEG) with the short- 
range electron-electron interaction w^ = w ee — w l J~ e . The 
analytic expression of e^ r (n) is known [4|| |47| and is re- 
called in Appendix IB II The short-range adiabatic LDA 
exchange kernel is given by the second-order derivative 
of the energy density with respect to the density, 



/^LDAl r > r ) 



d 2 t 



(n(r)) 



di 



5{r - r') 



(22) 



Just like its full-range LDA counterpart, the short-range 
exchange LDA kernel is thus strictly local in space. How- 
ever, this is here a less drastic approximation than for 
the full-range case. Indeed, by using the asymptotic 
expansion of the exact short-range spin-independent ex- 



change density functional for \x — > oo [l|, |43|, E & J\ 
— 7r/(4/x 2 ) J n(r) 2 dr+0 (l//i 4 ) , one can see that the exact 
adiabatic short-range exchange kernel has the following 
expansion in 1/fi, 



/r(r.r') 



2/z 2 V 



O 



If 



(23) 



i.e., in the limit of a very short-range electron-electron in- 
teraction, it also becomes strictly local. More than that, 
the short-range LDA kernel of Eq. (|22|) is exact for the 
leading term of Eq. (j2"5| . as shown in Appendix IB II 

The short-range LDA exchange kernel for a fixed value 
of the range-separation parameter fi = 0.4 is shown in 
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FIG. 1: Second-order derivatives of the full-range (full 

line) and short-range (/i = 0.4, dashed line) LDA 

exchange energy density with respect to the density n 

as functions of the Wigner-Seitz radius r s . 



Fig. [I] as a function of the Wigncr-Scitz radius r s = 
(3/ (47m)) and compared with the full-range LDA ex- 
change kernel. The LDA exchange kernel is always nega- 
tive, which is a consequence of the concavity of the LDA 
exchange energy density curve as a function of the den- 
sity n. For high enough densities such that r s <C l//i, 
the short-range LDA exchange kernel reduces to the full- 
range one (see Appcndix lB 1[) . For larger values of r s , the 
short-range LDA exchange kernel is reduced compared to 
the full-range one, and, in the low-density limit r s — > oo, 
it tends to the finite value of —ir/2/j, 2 while the full-range 
LDA exchange kernel diverges to — oo. 



B. Correlation kernel 

The short-range spin-dependent LDA correlation en- 
ergy functional is written as 



^ldaKH = / er(n(r),m(r))dr, (24) 
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FIG. 2: Second-order derivatives of the full-range (full 

line) and short-range (/i = 0.4, dashed line) LDA 

correlation energy densities with respect to the density 

n (top) and to the spin magnetization m (bottom) 

evaluated at in = as functions of the Wigner-Seitz 

radius r.. 



where e^, r (n,m) = ne c (n,m) — ne l £(n,m) is the com- 
plement short-range correlation energy density, obtained 
from the correlation energy per particle of the stan- 
dard homogeneous electron gas (HEG), e c (n, m), j49[ 
and the correlation energy per particle of the HEG with 
the long-range electron-electron interaction, ^(n, m), as 
parametrized fro m q uantum Monte Carlo calculations 
by Paziani et al. [50|. Its expression is recalled in Ap- 
pendix [B2l The singlet and triplet short-range adiabatic 
LDA correlation kernels are local functions given by the 
second-order derivatives of the energy density with re- 
spect to the density n and the spin magnetization to, 
respectively, 



7 c s Wr,r') 



/c,LDA( r i r ) 



dn 2 
3 2 ef (n(r),m(r)) 



J(r-r'), (25) 



S(t-t'). (26) 



dm 2 

For closed-shell systems, these kernels need to be 
evaluated only at zero spin magnetization, to = 0. 
Again, it can be argued that the strictly local form 
of the LDA correlation kernels of Eq. (|25j) and ([26]) 
is more appropriate for the short-range kernels than 
for the full-range ones. Using the asymptotic expan- 
sion of the exact short-range correlation functional for 
(j, -> oo 0, HI, E s c r [n,m] = tt/(2^ 2 ) / n 2 , c (r, r)dr + 
2v / 27r/(3/i 2 )/n2(r,r)dr + O (l//x 4 ), and the total and 
correlation on-top pair densities in the strong-interaction 
limit of the adiabatic connection A — > oo (or for fully 
spin-polarized systems n = \m\) [52l l53j. n 2 (r, r) — > 
and ri2,c(r, r) — > — ?i(r) 2 /2 + m(r) 2 /2, it is easy to show 
that the leading terms in the expansions of the exact adi- 
abatic short-range correlation kernels for /i — > oo, in the 
strong- interaction (or low-density) limit, are strictly local 



7f(r,r') 



2ft 



■6(r-r') + 



1 



3 £sr 



f c r (r,r') 



A^oo 2/X 2 



<5(r 



O 



(27) 



(28) 



The short-range LDA correlation kernels of Eqs. ([2"5"[) 
and ([26]) . using the parametrization of Ref. [50, are ex- 
act for these leading terms. 

The singlet and triplet short-range LDA correlation 
kernels are plotted in Fig. [51 and compared with the full- 
range LDA correlation kernels. The singlet LDA corre- 
lation kernel is always negative while the triplet LDA 
correlation kernel is always positive, reflecting the fact 
that the LDA correlation energy density is concave when 
plotted as a function of the density n and convex when 
plotted as a function of the spin magnetization to. As for 
the exchange kernels, the singlet and triplet short-range 
LDA correlation kernels reduce to the full-range kernels 
in the high-density limit r s — > (see Appendix IB 2 [I . In 
the low-density limit r s — > oo, they tend to the finite 
values of =F7r/2/i 2 , while the full-range kernels diverge to 
=Foo. 



IV. COMPUTATIONAL DETAILS 

The spin-adapted linear-response RSH scheme with 
the short-range LDA kernels has been implemented in 
a development version of the quantum chemistry pro- 
gram MOLPRO [55| for closed-shell systems. The im- 
plementation includes as special cases: standard TDKS 
with the LDA exchange-correlation functional, and time- 
dependent Hartree-Fock (TDHF). The implementation 
also includes the possibility to perform linear-response 
LC calculations (with the full-range LDA correlation 
functional). Each calculation is done in two steps: a 
self-consistent ground-state calculation is first performed 
with a chosen energy functional, and then a linear- 
response excited-state calculation is performed with a 
chosen kernel and using the previously calculated or- 
bit als. 

For compactness, "TD" will be dropped in the names 
of the methods and "LDA" will also be omitted in 
the names as it is the only density functional used 
here. Therefore, "KS" will denote a TDKS calculation 
using the LDA exchange-correlation functional, "HF" 
will stand for a TDHF calculation, "RSH" will de- 
note a linear-response RSH calculation using the short- 
range LDA exchange-correlation functional, and "LC" 
will stand for a linear-response LC calculation using the 
short-range LDA exchange functional and the full-range 
LDA correlation functional. We will call "RSH-TDA" a 
linear-response RSH calculation with the Tamm-Dancoff 
approximation. For all these methods, the same func- 
tional is used for the ground-state energy calculation and 
the linear-response calculation. In addition, we will call 
"RSH-fHx" a linear-response RSH calculation where only 
the Hartree-exchange part of the RSH kernel is used (no 
correlation kernel) but evaluated with regular RSH or- 
bitals (including the short-range correlation energy func- 
tional). 

Wc calculate vertical excitation energies and oscillator 
strengths of five small molecules, N 2 , CO, H 2 CO, C 2 H 4 , 
and QqH§, which have already been extensively studied 
theoretically [U, |H HHHI and experimentally [IIMI- 
In order to have unique, comparable references, equation- 
of-motion coupled-cluster singles doubles (EOM-CCSD) 
calculations were done in the same basis with the quan- 
tum chemistry program Gaussian 09 [64j . For each 
molecule, we report the first 14 excited states found with 
the EOM-CCSD method. Following Ref. HI we define 
the coefficient of the (spin-orbital) single excitation i — > a 
in the wave function of the excited state n to be 



1 



(A-B) 1/2 z r 



(29) 



but other choices for analyzing the eigenvectors are pos- 
sible, such as defining the weight of the single excitation 
i — > a to be w n ^ a = X 2 ia — Y 2 ia [65|. Each excited 
state was thus assigned by looking at its symmetry and 
the leading orbital contributions to the excitation. When 
several excited states of the same symmetry and the same 
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FIG. 3: Singlet (a) and triplet (b) excitation energies of N2 (in eV) with respect to the range-separation parameter 

\i (in bohr" 1 ) calculated by the linear-response RSH method with the short-range LDA exchange-correlation 

functional at the equilibrium geometry [54J and with the Sadlej+ basis set. 



leading orbital contributions were obtained, the assign- 
ment was done by increasing order of energy. Some as- 
signments for C2H4 and CgHg were difficult and are ex- 
plained in Tables |IV] and El The Sadlej basis sets [66|,|67| 
were developed to describe the polarizability of valence- 
like states. As the description of Rydbcrg states requires 
more flexibility, they were augmented with more diffuse 
functions to form the Sadlej 4- basis sets [56| that we use 
here. The molecules are taken in their experimental ge- 
ometries 0,111423. 

The C2H4-C2F4 dimer [H, [ill, [zl was studied in its 
cofacial configuration along the intermolecular distance 
coordinate R in the standard 6-31 G* basis set. A geome- 
try optimization was performed during the self-consistent 
ground-state calculation for each method. The CT exci- 
tation was identified by assigning the molecular orbitals 
involved in the excitations either to C2H4 or C2F4, using 
the visualization program MOLDEN [73| . 



V. RESULTS AND DISCUSSION 

Variation of the RSH excitation energies with 
the range-separation parameter 



> 



< 




0.2 0.3 0.4 
/i in bohr -1 



0.5 0.6 



FIG. 4: Mean absolute deviation (MAD) in eV of the 

first 14 excitation energies of the N2, CO, H2CO, C2H4 

and CgHg molecules calculated by the linear-response 

RSH method with the short-range LDA 

exchange-correlation functional with respect to the 

EOM-CCSD reference as a function of the 

range-separation parameter \x. 



When the range-separation parameter is zero, /1 = 0, 
the long-range HF exchange vanishes and the short-range 
exchange-correlation functional reduces to the usual full- 
range one, therefore the RSH method is equivalent to 
the standard KS method in this limit. With the LDA 
functional, linear-response KS gives good results for the 
low-lying valence excitation energies but underestimates 
the high-lying Rydberg excitation energies. This under- 
estimation is known to be due to the incorrect exponen- 
tial asymptotic decay of the LDA exchange potential [56| . 



When /1 increases, long-range HF exchange replaces LDA 
exchange and long-range LDA correlation is removed. In 
the limit /i — > +00, RSH becomes equivalent to a HF cal- 
culation, in which Rydberg excitation energies are usu- 
ally better described than in LDA but valence excitation 
energies can be poorly described, especially the triplet 
ones which can be strongly underestimated and can even 
become imaginary due to instabilities (A ± B in Eq. (TT4]) 
are no longer positive definite). 
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FIG. 5: Singlet (a) and triplet (b) excitation energies of N 2 , CO, H 2 CO, C2H4 and CeH 6 calculated by 

linear-response HF and KS (with the LDA functional), by the linear-response range-separated method RSH (with 

the short-range LDA functional and fj, = 0.4 bohr -1 ), as compared with the EOM-CCSD reference calculations. 



The variation of the first few singlet and triplet RSH 
excitation energies of N 2 with respect to the range- 
separation parameter /i is shown in Figure [3] The evo- 
lution of the excitation energies is similar for both spin 
states, however three different trends are seen for these 
excitations depending on their valence or Rydbcrg char- 
acter and their spatial symmetry The excitation ener- 
gies to the Rydberg excited states ( 1 S+, 1 n„, 1 E+, 3 S+, 
3 S+, 3 n u ) which are underestimated in KS show a sig- 
nificant increase with \x for /i > 0.1. This behavior is 
quite independent of the spin and spatial symmetry of 
the state. The valence excited states ( 1 n g , 1 S^, 1 A M , 
3 £+, 3 n ff , 3 A U , 3 E~, 3 n„) which arc correctly described 
in KS are less affected by the introduction of long-range 
HF exchange. However, we observe two opposite behav- 
iors depending on the orbital character of the excitation: 
all the valence n states (corresponding to a — > it orbital 
transitions) have excitation energies that slowly increase 
with n, while for valence £ and A states (correspond- 
ing to 7T — > 7T orbital transitions) the excitation energies 
decrease with fi. As a consequence, the ordering of the 
states changes significantly with /i. One should note that 
the variation of the excitation energies with /x have two 
causes: the variation of the orbital eigenvalues with /i 
in the ground-state calculation, and the variation of the 
kernel with [i in the linear-response calculation. Both 
effects can be significant. 

The choice of the range-separation parameter /i is im- 
portant. It has been proposed to adjust the value of \x 
for each system by imposing a self-consistent Koopmans' 
theorem condition [74|, UM ■ This approach is appealing 
but it has the disadvantage of being non size-consistent, 
so we prefer to use a fixed value of //, independent of 
the system. In Figure [51 the mean absolute deviation 
(MAD) of the RSH excitation energies with respect to 



the EOM-CCSD reference is plotted as a function of /j, 
for each molecule and for the total set. The global min- 
imum is obtained around \i w 0.3 — 0.4 bohr . In all 
the following, we use a fixed value of \i = 0.4, which 
is identical or similar to the value used in other range- 
separated TDDFT approaches [H HI HI . We note how- 
ever that the fact that the minimum of the MAD for 
CgHg is around /i = 0.2 shows that the optimal value of 
fi can substantially depend on the system. In particular, 
the presence of a triplet near-instability favors smaller 
values of \x. 



B. Accuracy of the RSH excitation energies and 
oscillator strengths 

The excitation energies and oscillator strengths for 
each method and each molecule are given in Tables [I]- 
fVl Mean absolute deviations with respect to the EOM- 
CCSD reference are also given for the valence, the Ry- 
dberg and all the excitation energies. We also report 
the position of the ionization threshold for each DFT 
method, as given by the opposite of the HOMO orbital 
energy. The excitation energies for all molecules are also 
plotted in Figure \5\ As expected, KS gives reasonably 
small errors for the valence excitation energies (MAD be- 
tween 0.36 and 0.72 eV) but deteriorates for the Rydbcrg 
ones (MAD between 0.49 and 1.83 eV) which are largely 
underestimated, as seen in Figure [5] As well known [56J, 
in KS with the LDA functional, the ionization energy is 
much too small, resulting in most of the Rydberg states 
and some of the valence states being in the continuum 
above the ionization threshold, and which are thus very 
much dependent on the basis set. This problem is absent 
in HF and range-separated approaches which correctly 



push up the ionization threshold. The HF excitation en- 
ergies arc usually larger than the reference ones except 
for the first triplet excitation energies which are much 
too small or even imaginary because of the HF triplet 
(near-)instability. Overall, HF gives relatively large total 
MADs (between 0.59 and 1.62 cV). 

The RSH excitation energies are in general interme- 
diate between KS and HF ones and in good agreement 
with the EOM-CCSD ones. The valence and Rydberg 
excitation energies are treated with a more uniform ac- 
curacy (MAD between 0.06 and 0.61 eV). However, the 
first triplet excitation energies are affected by the HF 
triplet near-instability and can be very underestimated. 
This effect is particularly important for the first triplet 
excitation energy of C2H4 and CeHg as shown in Ta- 
bles IIVI and [V] This underestimation is largely cured 
by the Tamm-Dancoff approximation, as shown by the 
RSH-TDA results. The quality of the other excitation 
energies is not deteriorated with this approximation, so 
that RSH-TDA gives overall smallest MADs than RSH. 
However, the oscillator strengths which were relatively 
good in RSH tend to be overestimated for excitations 
to valence states in RSH-TDA. This has been connected 
with the fact that the TDA oscillator strengths violate 
the Thomas-Reiche-Kuhn sum rule. The present RSH 
results give thus very much the same trends already ob- 
served with other types of range-separated TDDFT ap- 
proaches [H M M, tzi [z3 . 

The first singlet CT excitation energy in the C2H4- 
C2F4 dimer along the intermolecular distance coordinate 
R, for R between 5 and 10 A (i.e. between 9.45 and 
18.90 bohr), is given in Fig. [6] This excitation corre- 
sponds to an electron transfer from the HOMO of C2F4 
to the LUMO of C2H4. Therefore, its energy must be- 
have asymptotically as Ic 2 f 4 ~Ac 2 h 4 — 1/-R, where Ic 2 f a 
is the ionization potential of the tetrafluoroethylene and 
A-dHi is the electron affinity of ethylene. A fit of the 
form a + b/R was performed and the fitted parameters 
are shown in Fig. [5] The well-known deficiency of KS 
with the LDA functional to describe the —1/R depen- 
dence of such excitations is observed as it gives a param- 
eter b close to zero, while HF and RSH both give the 
expected correct asymptotic behavior in — 1/R thanks to 
the non-locality of their exchange kernel [7l| . 



C. Effect of the LDA correlation 

Tables ITHVI also report results obtained with the LC 
scheme using the short-range exchange LDA functional 
and the full-range LDA correlation functional. The com- 
parison with the RSH results allows one to see the global 
effect of long-range LDA correlation in the ground-state 
calculation and in the linear-response kernel. The RSH 
and LC excitation energies are globally quite close to each 
other, the largest difference being of 0.2 eV for the 3 n 
Rydberg state of the CO molecule. In most cases, the LC 
excitation energies are slightly larger than the RSH ones. 
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In comparison to the RSH scheme, the LC scheme gives 
slightly smaller MADs (by 0.01 to 0.08 eV) for valence 
excitation energies, but with the exception of CO it gives 
larger MADs (by 0.07 to 0.09 eV) for Rydberg excitation 
energies. The RSH and LC oscillator strengths are quite 
similar. This shows that long-range LDA correlation has 
a quite small effect for the systems and states considered 
here, and can be disregarded without much consequence. 
The first CT excitation energy in the C2H4-C2F4 dimer 
obtained with the LC scheme is also reported in Fig. [6] 
Not surprisingly, the RSH and LC curves have the same 
— 1/R behavior, which is given by the long-range HF ex- 
change kernel, and are essentially on-top on each other, 
showing that long-range LDA correlation has almost no 
effect on the HOMO and LUMO orbital energies. 

To investigate the effect of the short-range LDA cor- 
relation kernel, Tables HWl also report RSH-fHx results 
obtained with regular RSH orbitals but no correlation 
kernel at all. Removing the short-range LDA correla- 
tion kernel tends to yield larger singlet excitation ener- 
gies and smaller triplet excitation energies. This is not 
unexpected since the singlet LDA correlation kernel is 
negative and the triplet LDA correlation kernel is posi- 
tive, as shown in Figure [2] In comparison to the RSH re- 
sults, RSH-fHx gives quite similar singlet valence excita- 
tion energies and Rydberg excitation energies, but much 
lower triplet valence excitation energies (sometimes by 
as much as 0.5 eV), leading to significantly larger MADs 
for valence excitations. The short-range part of the LDA 
correlation kernel is thus important and cannot be ne- 
glected. 



VI. CONCLUSION 



Appendix A: Spin-adapted kernels 



We have studied a linear-response range-separated 
scheme, which combines a long-range HF exchange kernel 
with a short-range LDA exchange-correlation kernel, for 
calculating electronic excitation energies and oscillator 
strengths of molecular systems. It is a first-level approxi- 
mation before adding an explicit treatment of long-range 
correlation. It can also been seen as an alternative to 
the widely used linear-response LC scheme which com- 
bines a long-range HF exchange kernel with a short-range 
DFT exchange functional and a full-range DFT correla- 
tion functional. 

Tests on the N 2 , CO, H 2 CO, C 2 H 4 , and C 6 H 6 
molecules have shown that a reasonable value for the 
range-separation parameter is /i — 0.4 bohr -1 , which is 
consistent with what was previously reported in the liter- 
ature for other types of range-separated TDDFT meth- 
ods. Just like in the LC scheme, the introduction of 
long-range HF exchange in the present method corrects 
the well-known underestimation of high-lying Rydberg 
excitation energies of standard TDDFT using (semi)local 
density-functional approximations, but also leads to un- 
derestimated excitation energies to low-lying spin-triplet 
valence states. This latter problem is known to be asso- 
ciated with the presence of HF triplet near-instabilities 
and is largely cured by the Tamm-Dancoff approxima- 
tion which leads to a relatively uniform accuracy for all 
excitation energies, but possibly at the cost of deterio- 
rating the oscillator strengths. As expected, tests on the 
first CT excitation energy in the C2H4-C2F4 have shown 
that the present range-separated TDDFT method also 
correctly describe this kind of excitations. 

For the systems and states considered here the pres- 
ence of long-range LDA correlation in the ground-state 
calculation and in the linear-response kernel has a quite 
small effect, so that the present method gives results very 
similar to the ones given by the LC scheme. Long-range 
LDA correlation can therefore be disregarded. In con- 
trast, the short-range LDA correlation kernel is impor- 
tant for singlet — > triplet valence excitation energies and 
cannot be neglected. This work thus suggests that the 
present range-separated TDDFT scheme is a reasonable 
starting approximation for describing electronic excita- 
tion energies. The next step of this work is then to add 
to the present method an explicit frequency-dependent 
long-range correlation kernel derived from perturbation 
theory, e.g. in the spirit of Refs.[78l.[79l which would add 
the possibility of describing double excitations. 
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For spin-restricted closed-shell calculations, spin- 
singlet and spin-triplet excitations can be decoupled [42| - 
|44| (see also Refs. H| I80I ). The non-spin- flip part of the 
coupling matrix K of Eq. (p~2]l has the following spin block 
structure 



_ I K t,t K t4 



(Al) 



where the matrix blocks K^^ have elements of the form 
K^cr aa, ja' ha' with i, j, and a, b referring to occupied and 
virtual spatial orbitals, respectively. The matrix K can 
be brought to a block diagonal form by rotation in spin 
space 



K 



with a singlet component 



X K 

3 K 



K 



ICt.T + K^ + K; : -|- + K44 



and a triplet component 



5 K 



K 



T,T 



K t ,i - K u + K u 



(A2) 



(A3) 



(A4) 



This directly leads to the singlet and triplet RSH cou- 
pling matrices of Eqs. (fTS]) and P^|) . where the singlet 
and triplet short-range exchange-correlation kernels are 
defined as 



1 £ ST ( 



/lc( r l' r 2) — /zc.ti^ 1 " 1 ' 1 " 2 ) + /sc.t4,( r l' r 2) 



+ fZ,i t (r u r 2 ) + fZ, u (r u r 2 ) /4 



(A5) 



and 



3 .psr ( 



/zc( r l' r 2) = /xc,tf( ri ' r2 ) - /xc,n( r l' r 2) 

/^4t( r i' r 2) + /^ u (r 1 ,r 2 )l/4 



(A6) 



The different spin components of the kernel can be ex- 
pressed with the second-order functional derivatives of 
the corresponding energy functional E^[n,m] with re- 
spect to the density n and the spin magnetization density 
m, 



/xc,ft( r l' r 2) 

_ S 2 EZ[n,n 



5 2 E s x r c [n,m] 
r5nf(ri)£nt(r 2 ) 
, „ S 2 EZ[n,m] 



Jn(ri)(5n,(r 2 ) 



(5n(ri)(5m(r 2 ) 5m (ri) 6m (r 2 ) 

' (A7) 



and 



/xc,n( r l' r 2) = /xc,lt( r l> r 2) 



5 2 EZ[ 



n, m\ 



ra, m] S 2 E^ r c [n,m] 



(A8) 



i5n(ri)(5n(r 2 ) (5m(ri)(5m(r 2 ) 



10 



and 



S 2 E s x r c [n,m] 
n S 2 EZ[n,m] S' 2 EZ[n,m] 



/Jc,u( r l' r 2) = 

_ 5 2 EZ[n,m] 

6n(ri)8n(r2) <5n(ri)<5m(r2) 



Jm(ri)5m(r 2 )' 

' (A9) 



The mixed derivative with respect to n and m cancels out 
in Eqs. (|A5[) and (|A6[) and we finally obtain the singlet 
and triplet kernels of Eqs. |Q7|) and (TT5|) . 



Appendix B: Short-range LDA exchange-correlation 
functional 



1. Short-range LDA exchange 

The short-range spin-independent LDA exchange en- 
ergy density is a function of the density n (or equivalcntly 
of the Wigner-Seitz radius r s = (3/ (47m)) 1 / 3 ) and of the 
range-separation parameter /i, and writes 



n(e x 



(Bl) 



where e x is the full-range exchange energy per particle of 
the homogeneous electron gas, 



27a 2 
16 r ' 



(B2) 



with a = (4/(97r)) 1 / 3 , and e^T is the long-range exchange 
energy per particle of the homogeneous electron gas [4y, 



Jr 



9a" A 


[^UO 




r s 




(2A- 


4A 3 ) e -V(^ 2 ) . 


- 3A + 4A 3 



(B3) 



with A = a fir s /2. 



a. Large fi behavior 

In the limit of a very short-range interaction (/x — > 
+oo) or in the low-density limit (n — > or r s — > oo), 
i.e. A — > oo, the short-range exchange energy density e x r 
goes to zero with the following asymptotic expansion 



3 n 



9n 



16 r-3 fi 2 320 a 2 r 5 s ii A 



O 



/'■' 



(B4) 



and the corresponding short-range exchange kernel, i.e. 
the second-order derivative with respect to n, expands as 



1 sr 



d e^ 
dn 2 



2fi 2 6a 2 r 2 fi 



2 r 2 ,,4 



O 



1 



(B5) 



b. Small fi behavior 

In the limit of the Coulombic interaction (/i — > 0) or in 
the high-density limit (n — > +oo or r s — > 0), i.e. A — > 0, 
the short-range exchange energy density e x r reduces to 
the full-range exchange energy density e x = ne x with 
the following expansion 



fin 3ar s fi 2 n fi 4 
yfK 2tt 6tt 3 



O e 



-l/H' 



and the short-range exchange kernel behaves as 
d 2 ei r d 2 e x 



4 4 2 

7r a r s fi 



O 



dn 2 dn 2 
with the full-range exchange kernel 

8 2 e 



-l/ti- 



dn 2 



2 2 

— ir a r c . 



(B6) 



(B7) 



(B8) 



Taking the ratio of Eqs. (|B7[) and (|B8|) . it is seen the 
short-range exchange kernel reduces to the full-range one 
when 



a 2 fi 2 r 2 < 1, 



(B9) 



i.e. r s < 4.8 for /i = 0.4. 



2. Short-range LDA correlation 

The short-range spin-dependent LDA correlation en- 
ergy density is a function of the density n, the spin mag- 
netization density m (or equivalcntly of r s and Q = m/n) , 
and of the range-separation parameter (i, and writes 



n(eo-<), 



(B10) 



where e c is the full-range correlation energy per particle 
of the homogeneous electron gas [49| , and ej, r is the cor- 
relation energy per particle of a homogeneous electron 
gas with the long-range electron-electron interaction, as 
fitted by Paziani et al. [5(| on quantum Monte Carlo cal- 
culations with imposed exact limits, 



MCY 



>(C) 



+ ai(r s ,C)Ai 3 



+ a 2 (r s , Ofi 4 + a 3 (r s , ()fi 5 + a 4 (r s ,0fi 6 
+ a 5 (rs,C)^ 



(Bll) 



(l + b 2 (r s )f, 2 



where <fc>(0 = [(1 + C) 2/3 + fl - C) 2/3 ]/2 and the other 
functions are given in Ref. |5Q| . 

The derivatives of e s J with respect to n and m are 
easily expressed in terms of the derivatives of e s J with 
respect to r s and £. The first-order derivatives are 



de" 



r s de s n 



dn 

de s / 



3 dr 



def. 



(B12) 



9to d( 
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and the second-order derivatives are 



2 sr 



d 2 e 



dn 2 
d 2 e s J 



9n \ dr s 



dr 2 



(B13) 



dm 2 n dC, 2 



For spin-restricted closed-shell calculations, we just need 
to evaluate them at C = 0. 



a. Large \i behavior 

The leading terms of the asymptotic expansion for 
fi — > +00 of the short-range correlation energy density 
e s c r can be expressed with the on-top pair-density of the 
homogeneous electron gas [50|. In the low-density limit 
r s — > +00 (or the strong-interaction limit A — > +00 of 
the adiabatic connection), it simplifies to 



>+oo 



3(1-C 2 ) 
16 rf /j, 2 



O 



(B14) 



In this limit, the associated singlet and triplet short- 
range correlation kernels, i.e. the second-order deriva- 
tives of e^ r with respect to n and m evaluated at £ = 0, 
have the following expansions 



dn 2 



C— 0, r s — 7-+00 



■K ( 1 

v +0 b 



2 c sr 



8 2 e 



dm 2 



C=0, r s ^+oc 



7T 



o 



(B15) 



(B16) 



b. Small fi behavior 

In the limit /i — > 0, the short-range correlation energy 
density e* r reduces to the full-range correlation energy 
density e c =ne c with the following expansion [50[ 



e„ = e c 



2^ 



0{^). 



(B17) 



It can easily be shown that the singlet and triplet short- 
range correlation kernels, evaluated at C, = 0, approach 



the corresponding full-range kernels with the same lead- 
ing term 



2 c sr 



d 2 e 



dn 2 



dm 2 



C=o 



C=o 



dn 2 



4 4 2 

-ira r s fi 



0(^ 3 ), (B18) 



dm 2 



C=o 



C=o 



-Tra 4 r 4 s n 2 + 0{fi 3 ). (B19) 



In the high-density limit r s —J- 0, the expansion of the 
full-range correlation energy density has the form [49( 



e c = n[c (C)lnr s - ci(C) +0(r s lnr s )] . 



(B20) 



The expansion of the singlet full-range correlation kernels 
is 



d 2 e c 



dn 2 



Co(0) ir 2 a 3 r 3 s +0(r 4 s In r s ), (B21) 



C=o 



with co(0) = (1 — ln2)/7r 2 , and the expansion of 
the triplet full-range correlation kernels is found from 
the the correlation part of the spin stiffness a c (r s ) = 
(d 2 e c /d( 2 k=o 



8 2 e c 



dm 2 



C=o 



37r a r s a c (r s ) 



37H a 6 r 3 [A a In r s + C a + 0(r s In r s )] 



(B22) 

where A a = -1/(6tt 2 ) and C a = 0.0354744 [p. Com- 
paring Eqs. (|B18j) and (|B21|) . it is seen the singlet short- 
range correlation kernel reduces to the full-range one 
when 



7r a r s /j, 
1 - In 2 



<1, 



(B23) 



i.e. r s ^C 1.2 for /i = 0.4, and, comparing Eqs. (|B19p 
and (|B22|) . the triplet short-range correlation kernel re- 
duces to the full-range one when 



ar s ^ 



3ir{A a lnr s + C a ) 
i.e. r s < 2.4 for fi = 0.4. 



«1, 



(B24) 
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State 


Transition 


KS RSH 


RSH-TDA LC RSH-fHx 


HF 


EOM-CCSD 




Valence excitation energies (cV) 


o s + 


l7r« 


-> l7Tg 


7.87 7.19 


7.41 


7.31 6.65 


3.47 


7.72 


3 n 


3<x g 


-> lrr a 


7.54 7.84 


7.87 


7.88 7.67 


7.62 


8.16 


3 A U 


In* 


-* l7r s 


8.82 8.26 


8.39 


8.31 8.03 


5.86 


9.07 


X 


3<x g 


-> l7T 3 


9.05 9.43 


9.53 


9.43 9.47 


9.77 


9.55 


3 s~ 


lTTu 


-* 1t 9 


9.65 9.23 


9.26 


9.22 9.23 


7.94 


10.00 


l vz 


1-7T„ 


-> It? 


9.65 9.23 


9.26 


9.22 9.23 


7.94 


10.24 


X A U 


lir u 


-> l^g 


10.22 9.90 


9.84 


9.90 9.95 


8.78 


10.66 


3 n„ 


2er„ 


-» 1t s 


10.36 10.77 


10.81 


10.82 10.53 


11.28 


11.36 


Rydbcrg excitation energies (cV) 


>K 


3<Tg 


-*■ 4ct b 


10.28 11.78 


11.78 


11.94 11.73 


13.05 


11.74 


L K 


3cr, 


-> 4er g 


10.39 12.26 


12.29 


12.38 12.26 


13.98 


12.15 


3 s+ 


3er g 


-> 3cr u 


10.62 12.63 


12.63 


12.87 12.59 


14.16 


12.70 


3 n u 


3c g 


-> 2tt„ 


10.99 12.62 


12.62 


12.83 12.59 


14.56 


12.71 


l n u 


3c 9 


-> 2tt„ 


10.98 12.74 


12.74 


12.87 12.75 


13.21 


12.77 


x s+ 


3<X g 


-> 3cr„ 


10.62 12.76 


12.77 


12.89 12.78 


14.00 


12.82 








Ionization threshold: — 


ehomo (eV) 












10.38 15.34 


15.34 


15.76 15.34 


16.74 






MAD of excitation energies 


with res 


poet to EOM-CCSD (cV) 




Valence 






0.49 0.61 


0.55 


0.58 0.75 


1.82 


- 


Rydbcrg 






1.83 0.06 


0.07 


0.15 0.07 


1.35 


- 


Total 






1.06 0.38 


0.34 


0.40 0.46 


1.62 


- 








Oscillator 


strengths 


(XIO-*) 






x n„ 


3ct 9 


-> 2tt„ 


2.41 9.49 


9.54 


12.77 9.00 


8.42 


8.51 


x s+ 


3cr s 


-J- 3cr u 


1.06 21.11 


20.09 


27.59 19.94 


73.31 


17.36 



TABLE I: Excitation energies and oscillator strengths of N2 calculated by linear-response HF and KS (with the 

LDA functional), by the linear-response range-separated methods RSH, RSH-TDA, LC, and RSH-fHx (with the 

short-range LDA functional and ijl = 0.4 bohr -1 ), and by EOM-CCSD taken as reference, using the Sadlej+ basis 

set. Excitation energies above the ionization threshold are indicated in italics. 
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I State 



Transition 



KS 



RSH RSH-TDA LC RSH-fHx HF EOM-CCSD 









Valence 


excitatior 


energies 


(eV) 








3 n 


5ai(cr 


— > 2ei(7r*) 


5.95 


5.95 


6.02 


6.06 


5.62 


5.28 


6.45 


3 E + 


lei (7T 


-> 2ei(7r*) 


8.38 


8.15 


8.28 


8.22 


7.72 


6.33 


8.42 


1 n 


5ai(cr 


-> 2ei(7T*) 


8.18 


8.49 


8.67 


8.49 


8.55 


8.80 


8.76 


3 A 


lei(7r 


-> 2ei(7r*) 


9.16 


8.99 


9.07 


9.02 


8.80 


7.87 


9.39 


3 S~ 


lei(7r 


-> 2ei(7r*) 


9.84 


9.77 


9.79 


9.75 


9.77 


9.37 


9.97 


1 T,- 


lei (tt 


— > 2ei(ir*) 


9.84 


9.77 


9.79 


9.75 


9.77 


9.37 


10.19 


'A 


lei (tt 


— > 2ei(ir*) 


10.31 


10.31 


10.27 


10.29 


10.35 


9.96 


10.31 


3 n 


4ai(cr 


-> 2ei(7T*) 


11.40 


12.05 


12.07 


12.05 


11.91 


13.05 


12.49 






Rydbcrg 


excitation energies (cV) 








3 £+ 


5ai(<7 


— > Gai(a) 


9.55 


10.55 


10.56 


10.72 


10.46 


11.07 


10.60 


x £+ 


5ai(cr 


— > 6ai(o-) 


9.93 


11.32 


11.37 


11.36 


11.34 


12.23 


11.15 


3 £+ 


5ai(cr 


-> 7ai(cr) 


10.26 


11.34 


11.35 


11.51 


11.29 


12.40 


11.42 


*E+ 


5ai(cr 


-> 7oi(ct) 


10.47 


11.58 


11.59 


11.63 


11.60 


12.78 


11.64 


3 n 


5ai(cr 


-> 3ei(7r) 


10.39 


11.53 


11.54 


11.73 


11.49 


12.52 


11.66 


L n 


5ai(er 


-> 3ei(7r) 


10.48 


11.72 


11.72 


11.81 


11.73 


12.87 


11.84 



Ionization threshold: -(homo (°V) 



9.12 



13.83 13.83 14.27 13.83 15.11 



MAD of excitation energies with respect to the EOM-CCSD calculation (cV) 



Valence 


0.36 


0.31 


0.25 


0.29 


0.45 


0.89 


- 


Rydbcrg 


1.21 


0.10 


0.10 


0.09 


0.13 


0.92 


- 


Total 


0.73 


0.22 


0.19 


0.20 


0.31 


0.91 


- 







Oscillator strengths (xlO z ) 








'n 


5ai(er) -^ 2ei(7r*) 


8.69 8.64 12.96 8.73 


8.49 


8.55 


8.66 


*E+ 


5ai(cr) — y 6ai(er) 


1.84 4.26 4.44 3.70 


3.85 


10.58 


0.58 


*E+ 


5ai(er) — > 7ai(a-) 


12.53 13.73 14.67 15.86 


13.73 


9.39 


20.71 


1 n 


5ai(er) — > 3ei(7r) 


2.71 4.72 4.36 5.45 


4.58 


5.13 


4.94 
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TABLE II: Excitation energies and oscillator strengths of CO calculated by linear-response HF and KS (with the 

LDA functional), by the linear-response range-separated methods RSH, RSH-TDA, LC, and RSH-fHx (with the 

short-range LDA functional and [i = 0.4 bohr -1 ), and by EOM-CCSD taken as reference, using the Sadlej+ basis 

set. Excitation energies above the ionization threshold are indicated in italics. 
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I State 



Transition 



KS 



RSH RSH-TDA LC RSH-fHx HF EOM-CCSD 









Valence 


excitation 


energies 


(cV) 








3 A 2 


2b 2 (n 


-> 26i(tt*) 


3.06 


3.17 


3.19 


3.16 


3.08 


3.44 


3.56 


L A 2 


2b 2 (n 


-►2&i(7r*) 


3.67 


3.84 


3.85 


3.82 


3.86 


4.41 


4.03 


a A, 


\b\ (7T 


-» 2bi(0 


6.22 


5.65 


5.87 


5.74 


5.25 


2.15 


6.06 


a Br 


5ai (ct 


) -» 2bl(7T*) 


7.74 


8.11 


8.14 


8.11 


7.99 


8.19 


8.54 






Rydbcrg 


excitation energies 


(eV) 








»B 2 


26 2 (n 


— ► 6ai(er) 


5.84 


7.06 


7.07 


7.17 


7.01 


8.09 


6.83 


L B 2 


2b 2 (n 


— > 6ai((r) 


5.92 


7.26 


7.26 


7.30 


7.28 


8.55 


7.00 


6 B 2 


2b 2 (n 


— > 7ai(a) 


6.96 


7.91 


7.92 


7.99 


7.86 


8.98 


7.73 


*A 1 


2b 2 (n 


-> 36 3 (<r) 


6.73 


8.01 


8.01 


8.17 


7.96 


9.19 


7.87 


l B 2 


2b 2 (n 


— » 7ai((r) 


7.04 


8.15 


8.16 


8.17 


8.17 


9.39 


7.93 


L A t 


2b 2 (n 


-> 36 3 (<r) 


6.77 


8.18 


8.19 


8.27 


8.19 


9.28 


7.99 


l A 2 


2b 2 (n 


-> 3bi(7r) 


7.55 


8.58 


8.58 


8.67 


8.58 


10.04 


8.45 


6 A 2 


2b 2 (n 


-> 3bi(7r) 


7.58 


8.57 


8.57 


8.70 


8.56 


9.84 


8.47 


i B 2 


2b 2 (n 


— > Sai(a) 


7.97 


9.12 


9.14 


9.24 


9.06 


10.24 


8.97 


l B 2 


2b 2 (n 


— > 8ai((j) 


8.17 


9.42 


9.43 


9.49 


9.44 


10.84 


9.27 



Ionization threshold: — chomo (cV) 



6.30 



10.63 11.06 10.63 12.04 



MAD of excitation energies with respect to EOM-CCSD (eV) 



Valence 


0.46 


0.36 


0.28 


0.34 


0.50 


1.19 


- 


Rydbcrg 


1.00 


0.18 


0.18 


0.26 


0.16 


1.39 


- 


Total 


0.84 


0.23 


0.21 


0.29 
„-2% 


0.26 


1.34 


- 



l B 2 


2b 2 (n) -> 6ai(er) 


3.13 


1.77 


2.09 


1.88 


1.69 


2.99 


2.15 


l B 2 


2b 2 (n) -> 7ai(cr) 


2.05 


4.58 


5.54 


5.04 


4.47 


4.46 


4.12 


L A 1 


2b 2 (n) -> 36 2 (ct) 


4.34 


5.35 


6.10 


6.07 


5.18 


21.31 


5.70 


l B 2 


2b 2 (n) -> 801 (<r) 


4.27 


4.01 


4.84 


4.08 


3.88 


6.65 


4.22 
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TABLE III: Excitation energies and oscillator strengths of H2CO calculated by linear-response HF and KS (with the 

LDA functional), by the linear-response range-separated methods RSH, RSH-TDA, LC, and RSH-fHx (with the 

short-range LDA functional and /i = 0.4 bohr -1 ), and by EOM-CCSD taken as reference, using the Sadlej+ basis 

set. The molecule is oriented in the yz plane along the z axis. Excitation energies above the ionization threshold arc 

indicated in italics. 
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I State 



Transition 



KS 



RSH RSH-TDA LC RSH-fHx HF EOM-CCSD 









Valence excitation 


energies 


(eV) 








»B lu 


li>3u(w 


-> lf>2g(0 


4.63 


3.78 


4.13 


3.92 


3.31 


0.16 


4.41 


L B Ul 


lb 3u (7T 


-> lf>2,(0 


749 


7.60 


8.03 


7.60 


7.65 


7.35 


8.00 


Big 


U 3g (a) 


-> U>2 B (lT*) 


7.18 


8.03 


8.04 


8.15 


8.02 


8.36 


8.21 


1 B lg 


lb 3g (a 


-> 162 9 (tt*) (a) 


7.47 


8.15 


8.16 


8.23 


8.16 


9.36 


8.58 




Rydbcrg 


excitation 


energies 


(eV) 








B3 u 


lb 3u (TV 


— ► 4ai g (ir) 


6.60 


7.32 


7.33 


7.44 


7.29 


6.87 


7.16 


l B Xu 


lb 3u (7T 


-> 4ai g (cr) 


6.67 


7.49 


7.49 


7.55 


7.50 


7.13 


7.30 


-Big 


lb 3u (TV 


-> 2h 2u (<r) 


6.98 


7.51 


7.53 


7.50 


7.43 


7.63 


7.91 


B 2g 


l&3u(w 


-^ 3hi„(<r) 


6.98 


8.16 


8.16 


8.26 


8.14 


7.75 


7.93 


'Big 


lb 3u (ir 


-> 2h 2u (cr) <b) 


7.20 


8.04 


8.05 


8.02 


8.05 


7.74 


7.97 


B^g 


lb 3u (7r 


-> 3hi u (<r) 


7.03 


8.27 


8.27 


8.36 


8.27 


7.91 


8.01 


3 A 


lb 3u (7T 


-> 26 3u (7r) 


8.08 


8.55 


8.55 


8.77 


8.49 


7.97 


8.48 


l A q 


lb 3 u(TT 


-> 26 3u (7r) 


8.32 


8.95 


8.98 


9.04 


8.97 


8.57 


8.78 


d B 3 „ 


lb 3u (7T 


— > 5a-i g (cr) 


8.09 


9.08 


9.09 


9.18 


9.03 


8.71 


9.00 


l B Su 


1&3uO 


— > 5a-i g (cr) 


8.10 


9.22 


9.23 


9.26 


9.23 


8.92 


9.07 



Ionization threshold: — chomo (°V) 



[ 6.89 10.61 10.61 11.07 10.61 10.23" 



MAD of excitation energies with respect to EOM-CCSD (cV) 



Valence 


0.72 


0.41 


0.23 


0.33 


0.52 


1.46 


- 


Rydbcrg 


0.76 


0.18 


0.18 


0.26 


0.18 


0.24 


- 


Total 


0.74 


0.24 


0.20 


0.28 


0.27 


0.59 


- 










-a-i 









l fli« 


lb 3u 


(it) -> 16 2g (0 


30.64 


35.42 


53.86 


35.77 


35.10 


39.99 


36.29 


l B 3u 


lb 3u 


(n) -> 4ai g (cr) 


6.30 


7.64 


8.02 


8.14 


7.38 


9.08 


8.16 


L B 3u 


lb 3u 


(71) -> 5ai g (cr) 


0.02 


1.26 


1.30 


1.56 


1.14 


0.63 


0.61 



TABLE IV: Excitation energies and oscillator strengths of C2H4 calculated by linear-response HF and KS (with the 

LDA functional), by the linear-response range-separated methods RSH, RSH-TDA, LC, and RSH-fHx (with the 

short-range LDA functional and [i = 0.4 bohr -1 ), and by EOM-CCSD taken as reference, using the Sadlej+ basis 

set. The molecule is oriented in the yz plane along the z axis. Excitation energies above the ionization threshold are 

indicated in italics. ( a ' and iyV) : These two excitations mix heavily in LDA [82, |83[ and the leading orbital 

contribution to the excitation changes with the range-separation parameter. Adiabatic curves respect to /z were 

followed to do the assignment, with state (b) defined as the lowest 1 B\ g state with orbital transitions 

l&3g (c) — > 1629(1"*) and l&3 U (7r) — > 2&2m(c), and state (a) defined as the second lowest one. 
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State Transition 



KS 



RSH RSH-TDA LC RSH-fHx HF EOM-CCSD 









Valence 


excitation 


energies 


(cV) 








*B lu 


le lg (7T 


— > le 2 u(7r*) 


4.35 


3.37 


3.94 


3.49 


2.88 


- 


3.96 


3 E lu 


le lg (n 


-> le 2 „(7r*) 


4.69 


4.81 


4.83 


4.84 


4.69 


4.68 


4.90 


1 B 2 „ 


le lg (TT 


— > le2 u (7T*) 


5.20 


5.45 


5.55 


5.45 


5.47 


5.78 


5.15 


3 B 2u 


lei 9 (7r 


— > le2u(7T*) 


4.94 


5.02 


5.19 


5.05 


4.95 


5.02 


5.78 


1 B lu 


le lg (7T 


— > le2 u (7T*) 


5.97 


6.25 


6.44 


6.24 


6.29 


5.84 


6.52 


1 E lu 


lei g (7T 


— > le 2 u(7T*) 


6. SO 


7.14 


7.64 


7.13 


7.16 


7.34 


7.30 






Rydbcrg 


excitation 


energies 


(eV) 








E\g 


lei g (7r 


-> 4a l9 (<r) 


6.01 


7.00 


7.00 


7.10 


7.00 


6.46 


6.40 


E\g 


le lg (7T 


-> 4ai 9 (<r) 


6.03 


7.08 


7.08 


7.14 


7.09 


6.59 


6.46 


3 A 2u 


le lg (7T 


— > 4ei u (a) 


6.52 


7.43 


7.43 


7.56 


7.43 


6.87 


6.92 


'A 2u 


le lg (7T 


— > 4ei u ((r) 


6.5^ 


7.53 


7.53 


7.63 


7.53 


7.01 


7.00 


3 E 2u 


le lg (7T 


— > 4ei u (a) 


6.54 


7.68 


7.68 


7.83 


7.68 


7.17 


7.06 


1 E 2 „ 


le lg (7T 


— > 4ei u ((r) 


6.55 


7.71 


7.71 


7.84 


7.71 


7.21 


7.08 


'A lu 


le lg (7T 


— > 4ei u (a) 


6.59 


7.90 


7.90 


8.05 


7.90 


7.43 


7.18 


3 A lu 


lei 9 (7r 


— ► 4ei u ((r) 


6.59 


7.90 


7.90 


8.05 


7.90 


7.43 


7.19 



Ionization threshold: — chomo (cV) 



(i.oO 



9.72 



9.72 



10.18 



.72 



9.15 



MAD of excitation energies with respect to EOM-CCSD (eV) 



Valence 


0.39 


0.32 


0.22 


0.29 


0.43 


0.93 


Rydbcrg 


0.49 


0.61 


0.62 


0.74 


0.62 


0.12 


Total 


0.45 


0.49 


0.45 


0.55 


0.54 


0.47 






lei 9 (7r) -> le 2 „( 
leigjw) -> 4ei„( 



a) 



Oscillator strengths (xlO ) 



55.78 
2.11 



62.74 
7.10 



94.95 
7.55 



63.00 
8.27 



62.42 71.49 
6.87 7.69 



66.41 
7.04 



TABLE V: Excitation energies and oscillator strengths of C§Rq calculated by linear-response HF and KS (with the 

LDA functional), by the linear-response range-separated methods RSH, RSH-TDA, LC, and RSH-fHx (with the 

short-range LDA functional and /x = 0.4 bohr -1 ), and by EOM-CCSD taken as reference, using the Sadlcj+ basis 

set. Excitation energies above the ionization threshold are indicated in italics. Except for the E lg states whose 

assignments were trivial, all the other states correspond to the orbital transitions e\ g — > ei u and e\ g — > e\ u , which 

lead to B\ u E\ u B2 U and A\ u Ei u ® ^b« manifolds, respectively, in the 2?6/i symmetry point-group. Since the 

calculations were performed in the T>2h subgroup, some symmetry information for these states was lost but the 

assignment could be done using the degeneracy of the states. 



